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I. INTRODUCTION 

Models of gauge-mediated supersymmetry (SUSY) 
breaking (GMSB) pQ are attractive for several reasons. 
Theoretically they solve the "naturalness problem" [2 j 
and provide a low mass (warm) dark matter candi- 
date [3]. From an experimental standpoint they pro- 
vide a natural explanation for the observation of an 
eej-ffrr |H H] candidate event by the CDF experiment 
during Run I at the Fermilab Tevatron. In particular, 
the photon (7) and missing transverse energy (tfhr) can 
be produced by the decay of the lightest neutralino (x?) 
into a photon and a weakly interacting, stable gravitino 
(G). While much attention has been given to prompt 

X? 



tinue in Section VI with a description of the optimiza- 



7G decays, versions of the model that take into 



account cosmological constraints favor a G with keV/c 2 
mass and a Xi with a lifetime that is on the order of 
nanoseconds or more [6j. 

Here we describe in detail [7] the first search for heavy, 
long-lived neutralinos using photon timing at a hadron 
collider in the 7 + jet + frp final state where we require 
at least one jet and at least one photon. The data com- 
prise an integrated luminosity of 570±34 pb _1 of pp col- 
lisions at yfs = 1.96 TeV from the Tevatron collected 
with the CDF II detector [5]. Previous searches for 
sub-nanosecond [TU] and nanosecond-lifetime [TU] Xi 
— > 7G decays using non-timing techniques have yielded 
null results. The present results extend the sensitivity to 
larger \\ lifetimes and masses. 

The structure of this paper is as follows: the remain- 
der of this section provides a more detailed motivation for 
the search and describes the CDF detector, in particular 
the recently installed timing system on the electromag- 
netic calorimeters (the "EMTiming" system) that is used 
to measure the time of arrival of photons. Section [TtJ 
describes how photons from heavy, long-lived particles 
would interact with the detector and how the standard 
identification criteria for prompt photons are modified 
to keep the identification efficiency high for delayed pho- 
tons. The section further describes the photon timing 



measurement. We describe the data sample in Section III 
and discuss the event pre-selection criteria. Section IV 
describes the various background sources as well as the 
methods of estimating the rate at which they populate 
the signal region. After a description and estimation of 
the acceptance for GMSB events in Section W\ we con- 
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tion procedure and the expected sensitivity. The data 
are studied in Section IVIII and limits are set on GMSB 
with a model-independent discussion of the sensitivity. 
Section IVIIII concludes with the final results and a dis- 
cussion of the future prospects for a similar analysis with 
more data. 



A. Theory and Phenomenology 

Many minimal GMSB models are well specified with a 
small number of free parameters. The elcctroweak sym- 
metry breaking mechanism originates in a "hidden sec- 
tor" (not further specified in the model) and is medi- 
ated to the visible scalars and fermions by messenger 
fields; for more details see pQ and references therein. 
The free parameters of the minimal GMSB model are 
as follows: the messenger mass scale, M m ; the number 
of messenger fields, iV m ; a parameter A that determines 
the gaugino and scalar masses; the ratio of the neutral 
Higgs vacuum expectation values, tan(/3); and the sign 
of the higgsino mass parameter, sign(/i). For models 
with N m = 1 and low tan(/3) < 30 the weakly inter- 
acting G is the lightest supersymmetric particle (LSP) 
and the next-to-lightest supersymmetric particle (NLSP) 
is the lightest neutralino Xi- F° r models with N m > 1 
or tan(/3) > 30, the NLSP is a slepton (mostly f x ) [IT] . 
As there are many GMSB parameter combinations that 
match this phenomenology, representative "model lines" 
have been identified that allow a good specification of 
the model with only one free parameter that sets the 
particle masses. This analysis follows line 8 of the Snow- 
mass Points and Slopes (SPS 8) proposal [12] and as- 
sumes M m = 2A, tan(/3) = 15, sgn(/x) = 1, N m = 1, and 
R-parity conservation. In this model the Xi decays via 
Xi — * 7G with a branching ratio of ^100% but leaves the 
Xi mass and lifetime as free parameters. 

Non-minimal GMSB models with a non-zero Xi life - 
time and a ~ 1-1. 5 keV/c 2 mass G are favored as they are 
consistent with current astronomical observations and 
models of the early universe that take inflation into ac- 
count Q2]. If the G's are too light (<1 keV/c 2 ), they 
can destroy the nuclei produced during big bang nucle- 
osynthesis, leading to a cosmic microwave background 
that is different from observations [BJ. If they are too 
heavy (>1 keV/c 2 ), while they are a warm dark matter 
candidate [3] and consistent with models of galaxy struc- 
ture formation, their density can cause the universe to 
overdose. To include the proper GMSB messenger par- 
ticle decays and lifetimes, an additional SUSY breaking 
scale is included and provides an additional parameter 
in the model that relates the xi lifetime with the G and 
the Xi masses. In this formulation |1 our parameter 
choices, SPS 8, favor a lifetime of several nanoseconds 
for the 100 GeV/c 2 x? mass range, just above current 
exclusions [5J [TO] . 
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(b) 

FIG. 1: Feynman diagrams of the dominant tree-production 
processes at the Fermilab Tevatron for the SPS 8 GMSB 
model line. The taus and second photons, if available, can be 
identified as jets in the detector. Note that only one choice 
for the charge is shown. 

In pp collisions supersymmetric particles are pair- 
produced due to i?-parity conservation. We probe a 
range of A not already excluded at 95% confidence level 
(C.L.) in previous collider experiments 0HO] where the 
squarks and gluinos have masses of ^600-800 GeV/c 2 
and the sleptons and gauginos have masses of ~100- 
300 GeV/c 2 . At the Tevatron, with y/s = 1.96 TeV, 
squarks and gluinos are too heavy to have significant 
production cross sections, hence gaugino pair-production 
dominates pQ. Individually, X2X1 an d X1X1 production, 
as shown in Fig.[IJ contribute 45% and 25%, respectively, 
of the total GMSB production cross section (<T pro( j). The 
rest of the production is mostly slepton pairs. We note 
that (Tprod is independent of the \i lifetime. 

This analysis focuses on the 7 + Jfir final state which is 
expected to be more sensitive to the favored nanosecond 
lifetime scenario [2]. To identify GMSB events, we use 
the CDF II detector. As shown in Fig. [T] each gaugino 
decays (promptly) to a Xi i n association with taus whose 
decays can be identified as jets [15] . Whether the x? 
— > 7G decay occurs either inside or outside the detector 
volume depends on the Xi decay length (and the detec- 
tor size). The Xi' s and/or the G's leaving the detector 
give rise to frr since they are weakly interacting particles 



(the neutrinos in the event also affect the fk"). Depend- 
ing on whether one or two Xi s decay inside the detector, 
the event has the signature of high energy 77 + frp or 
7 + fr, often with one or more additional particles from 
the heavier sparticle decays. These are identifiable as an 
additional jet(s) in the detector. We do not require the 
explicit identification of a tau. This has the advantage 
of reducing the model dependence of our results, making 
them applicable to other possible gaugino decay mod- 
els. A study to see if there is additional sensitivity from 
adding r identification to the analysis is in progress. 

The arrival time of photons at the detector allows for 
a good separation between nanosecond- lifetime Xi' s an d 
promptly produced standard model (SM) photons as well 
as non-collision backgrounds. Figure |2ja) illustrates a Xi 
— > 7G decay in the CDF detector after a macroscopic 
decay length. A suitable timing separation variable is 

, ^ I <JC J- tJC ^ I / 1 \ 

^corr = \tf —ti) , (1) 

where tt — ti is the time between the collision ti and 
the arrival time tf of the photon at the calorimeter, and 
\xf — Xi\ is the distance between the position where the 
photon hits the detector and the collision point. Here, 
t colr is the photon arrival time corrected for the collision 
time and the timc-of- flight. Prompt photons will produce 
£ corr = while photons from long-lived particles will ap- 
pear "delayed" (i C orr > 0), ignoring resolution effects. 
Figure [2jb) shows the simulated distribution of t corr for 
a GMSB signal, prompt photons, and non-collision back- 
grounds in the detector. 

B. Overview of the Search 

This search selects photons with a delayed arrival time 
from a sample of events with a high transverse energy 
(E-p) isolated photon, large f!r, and a high- Erp jet to 
identify gaugino cascade decays. The background to 
this search can be separated into two types of sources: 
collision and non-collision backgrounds. Collision back- 
grounds come from SM production, such as strong inter- 
action (QCD) and electroweak processes. Non-collision 
backgrounds come from photon candidates that are ei- 
ther emitted by cosmic ray muons as they traverse the 
detector or are from beam related backgrounds that pro- 
duce an energy deposit in the calorimeter that is recon- 
structed as a photon. 

The search was performed as a blind analysis, pick- 
ing the final selection criteria based on the signal and 
background expectations alone. The background rates in 
the signal region are estimated using t colr control regions 
from the same 7 + jet + tf}j< data sample and comparing 
to the distribution shapes of the various backgrounds. 
A Monte Carlo (MC) simulation is used to model the 
GMSB event dynamics and timing in the detector and to 
estimate the signal expectations. Combining these back- 
grounds and signal event estimates permits a calculation 
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FIG. 2: (a) The schematic of a long-lived \i decaying into a 
G and a photon inside the detector. While the G leaves un- 
detected the photon travels to the detector wall and deposits 
energy in the detector. A prompt photon would travel di- 
rectly from the collision point to the detector walls. Relative 
to the expected arrival time, the photon from the x? would 
appear "delayed." (b) The i corr distribution for a simulated 
GMSB signal at an example point of m~o = 100 GeV/c 2 and 
t~o — 5 ns as well as for standard model and non-collision 
backgrounds. 



C. The CDF II Detector and the EMTiming 

System 

The CDF II detector is a general-purpose magnetic 
spectrometer, whose detailed description can be found 
in [8] and references therein. The salient components are 
summarized here. The magnetic spectrometer consists 
of tracking devices inside a 3-m diameter, 5-m long su- 
perconducting solenoid magnet that operates at 1.4 T. 
A set of silicon microstrip detectors (SVX) and a 3.1-m 
long drift chamber (COT) with 96 layers of sense wires 
measure the position [Si) and time (U) of the pp interac- 
tion and the momenta of charged particles. Muons from 
the collision or cosmic rays are identified by a system of 
drift chambers situated outside the calorimeters in the 
region with pseudorapidity \r]\ < 1.1. The calorimeter 
consists of projective towers (8<f> = 15° and Srj sw 0.1) 
with electromagnetic and hadronic compartments and is 
divided into a central barrel that surrounds the solenoid 
coil (| 77 1 < 1.1) and a pair of end-plugs that cover the 
region 1.1 < \rj\ < 3.6. Both calorimeters are used to 
identify and measure the energy and position of photons, 
electrons, jets, and fa. 

The electromagnetic calorimeters were recently instru- 
mented with a new system, the EMTiming system (com- 
pleted in Fall 2004), which is described in detail in [TB] 
and references therein. The following features are of par- 
ticular relevance for the present analysis. The system 
measures the arrival time of electrons and photons in 
each tower with \i]\ < 2.1 using the electronic signal 
from the EM shower in the calorimeter. In the region 
1 77 1 < 1.1, used in this analysis, photomultiplier tubes 
(PMTs) on opposite azimuthal sides of the calorime- 
ter tower convert the scintillation light generated by the 
shower into an analog electric signal. The energy mea- 
surement integrates the charge over a 132 ns timing win- 
dow around the collision time from ~20 ns before the 
collision until ^110 ns afterwards. New electronics in- 
ductively branches off ~15% of the energy of the anode 
signal and sends it to a discriminator. If the signal for a 
tower is above 2 mV (~3-4 GeV energy deposit), a digi- 
tal pulse is sent to a time-to-digital converter (TDC) that 
records the photon arrival time and is read out for each 
event by the data-acquisition system. The resolution of 
the time of arrival measurement is 0.50 ± 0.01 ns for the 
photon energies used in this analysis. 



II. PHOTON IDENTIFICATION AND TIMING 



of the most sensitive combination of event requirements. 
We note that the jet requirement helps make this search 
sensitive to any model that produces a large mass particle 
decaying to a similar final state. 



The CDF detector has been used for the identification 
(ID) of high-energy photons for many years, and a stan- 
dardized set of ID criteria (cuts) for the region \rj\ < 1.0 
is now well established. Each cut is designed to sepa- 
rate real, promptly produced photons from photons from 
7t° — > 77 decays, hadronic jets, electrons, and other back- 
grounds, see [3 |HJ [T7] for more details and the Appendix 
for a description of the ID variables. 
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^/CDF Calorimeter 




(a) 

CDF Calorimeter 




(b) 

FIG. 3: The definitions of the a and /3 incident angles using 
schematic diagrams of a long-lived Xi decaying to a photon 
and a G in the CDF detector. The angles a and f3 are the 
projections of the incident angle ip at the front face of the 
calorimeter in the (r,z)- and the (r,0)-plane, respectively. 



Unlike photons from SM processes, "delayed" photons 
from long-lived Xi ' s are n °t expected to hit the calorime- 
ter coming directly from the collision point [2]. As 
shown in Fig. |2^a), Xi' s with a long lifetime and small 
boost can produce a photon from Xi ~ * lG with a large 
path length from the collision position to the calorime- 
ter (large i C orr)- We define the photon incident angle at 
the face of the EM calorimeter, ip, as the angle between 
the momentum vector of the photon from the Xi an d the 
vector to the center of the detector. For convenience we 
consider the ip projection onto the (r, z)-plane and label 
it a, and the ip projection onto the (r, 0)-plane and label 
it /3; see Fig. [3j This distinction is made as the photon 
ID variable efficiencies vary differently between a and /3. 

Figure [4] compares the ip distribution for prompt, SM- 
like photons and photons from long-lived x?' s - Each 
are simulated as the decay product of a Xi with 
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FIG. 4: The distribution of the total incident angle tp at the 
front face of the calorimeter for simulated photons from x5's 
with m~o = 110 GeV/c 2 . "Prompt" photons from xV s with 

a lifetime of ns (solid) are compared to photons from Xi's 
with a lifetime 10 ns (dashed). The dotted histogram shows 
the distribution for a lifetime of 10 ns for photons with 2 < 
£ C orr < 10 ns and shows that, as expected, delayed photons 
can have a significant incident angle. 



7Ti~o = 110 GeV/c 2 using the PYTHIA MC generator [15] . 
The distributions of promptly produced photons |19j 
have a maximum at ip—0° and extend to ~18° in a while 
f3 is always <C 1° as the beam has negligible extent in the 
x-y plane. The most common angle ip for a simulated 
neutralino sample with r^o = 10 ns is ^10° and extends 
out to maximum angles of ~60° and ~40° in a and (3 
respectively. For this sample, the majority of photons 
arrive at angles between and 40° total incident angle. 
The mean of the distribution rises as a function of t^o 
but becomes largely independent of m^o and r^o in the 
range 10 < t^o < 35 ns. Also shown is the distribution 
for delayed photons, selected with 2 < t corr < 10 ns, sim- 
ilar to a typical final analysis requirement. The delayed 
photon requirement shifts the maximum of the distribu- 
tion of ip from ~10° to ^25°. As the incident angles of 
photons from long-lived particles are much larger than 
for prompt photons, the standard selection criteria are 
re-examined and modified where necessary. 

To verify that we can robustly and efficiently identify 
photons from heavy, long-lived particles, we examine the 
efficiencies of the photon ID variables as a function of a 
and j3 separately. As we will see the standard photon 
identification requirements are slightly modified for this 
search; each is listed in Table [I] To study photon showers 
at a wide variety of angles in the calorimeter, we create 
a number of data and MC samples of photons and elec- 
trons. An electron shower in the calorimeter is very sim- 
ilar to that from a photon, but electrons can be selected 
with high purity. We create two samples of W — > ev 
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Er > 30 GeV and \r)\ < 1.0 

Fiducial: not near the boundary, in <j> or z, of a 

calorimeter tower 
£ H ad/£ E M < 0.125 

Energy in a AR — 0.4 cone around the photon 
excluding the photon energy: 
E lso < 2.0 GeV + 0.02 • (Er - 20 GeV) 

No tracks pointing at the cluster or 

one with p T < 1.0 GeV/c + 0.005 ■ Er 

Ep T of tracks in a 0.4 cone < 2.0 GeV/c + 0.005 ■ Er 



Electron Requirements 



E 



2 na clustcr 



< 2.4 GeV + 0.01 • Er 



Ap = l-gPMTl-gpMT2l < g 



TABLE I: The photon identification and isolation selection 
requirements. These are the standard requirements with the 
Xces < 20 requirement removed. These variables are de- 
scribed in more detail in [5] [T7] and the Appendix. 



events, one from data, and the other simulated using the 
PYTHIA MC generator and the standard, GEANT based, 
CDF detector simulation [5D]. Each must pass the re- 
quirements listed in Table [TTJ Similarly, two samples of 
MC photons are generated using xl ~ ¥ lG decays with 
m^o = 110 GeV/c 2 and r^o = ns and r^o = 10 ns respec- 
tively to cover the region < ip < 60°. The highest Ep 
photon in the event is required to be the decay product 
of a Xi an d to pass the Er, 77, and fiducial requirements 
listed in Table III 

Figure [5] compares the distributions of the photon ID 



variables for the r T 



and t^o 



10 ns samples. A 



visual comparison snows that the differences are, on av- 
erage, very small. Figure [6] shows the efficiency of the 
final set of ID requirements for MC photons and elec- 
trons from data as a function of incident angles a and (3 
(taking the photon position from the measured center of 
the calorimeter energy cluster). The efficiencies are very 
similar and constant except at large values of (3 where 
the efficiency drops, which is where real collision data 
are not available. The drop in efficiency at large (3 is 
due to the photon shower in the calorimeter traversing 
into the neighboring tower in (p. Because the photons 
are identified and measured as clusters in the calorime- 
ter [17], this decreases the cluster-energy sum while in- 
creasing the isolation energy. Therefore, the photon ap- 
pears non-isolated and the isolation efficiency falls from 
-98% at (3 = 0° to -90% at (3 = 50°. This is not a 
problem for large a as energy leakage into the neighbor- 
ing tower in 77 is included in the energy sum. The total 
photon identification efficiency as a function of ip in this 
regime falls from —93% to —80%. However, since in our 
ip region the fraction of events with large j3 is small (see 
Fig.Wl, even at large r^o, the ID criteria are only —1.5% 



less el 



icicnt for photons for the t^o — 10 ns sample than 



for the prompt sample. Thus, the majority of the stan- 
dard requirements are not changed for the search. The 



Er > 30 GeV and H < 1.0 

Fiducial: not near the boundary, in <j> or z, of a 

calorimeter tower 
0.9 < E/p < 1.1 or p T > 50 GeV/c 
Track traverses >3 stereo and >3 axial COT 

superlayers with 5 hits each 
Additional requirements to reject electrons from 7 — > ee 
Global Event Requirements 

Er >30 GeV 

Exactly 1 vertex with iVtrks > 4 and \z\ < 60 cm 
Transverse mass of the electron and ftrr: 
50 < m T < 120 GeV/c 2 



TABLE II: The requirements used to select electrons from 
W — > ev events to validate the ID efficiency of simulated pho- 
tons. These are topological and global event cuts in combi- 
nation with loose calorimetry but tight track quality require- 
ments. This produces a sample that contains electrons with 
high purity but has a low bias for calculating the efficiency 
of photon ID requirements vs. incident angle. The vertex 
reconstruction algorithm is described in Secti on |II A| and uses 
tracks passing the requirements listed in Table [TTl| These vari- 
ables are summarized in the Appendix and described in more 
detail in [H]. 



efficiency variation as a function of angle is taken into ac- 
count by using the detector simulation for the efficiencies 
and assigning a 5% systematic uncertainty to the overall 
photon ID efficiency measurement. 

The comparison of the photon shower-maximum pro- 
file to test-beam expectations, Xces 03' ^ s removed 
from the photon identification requirements because it 
becomes inefficient at large angles. The shower for a pho- 
ton that hits the shower-maximum detector (CES) at a 
large value of a ((3) angle would spread out and have 
a larger-than-expected RMS in the z ((f)) direction due 
to the projection. A GEANT simulation [20 shows the 
efficiency of the Xces requirement is constant at small 
angles, but then falls off rapidly at large angles. Thus, 
the photon Xces requirement is removed. 

A second change to the standard photon ID is to add 
a requirement to remove high-energy photon candidates 
that are caused by a high- voltage breakdown ("spike") 
between the PMT photocathode and the surrounding 
material. Such an occurrence can produce false photon 
candidates that are uncorrelated with the collision and 
appear delayed in time. Spikes are identified by the asym- 
metry of the two energy measurements of the PMTs of a 
tower: 



\E 



PMTl 



E 



PMT2 I 



E 



PMTl 



E 



PMT2 



(2) 



where -Epmti an d -EpMT2 are the two PMT energies. Fig- 
ure[7]compares photon candidates from both real photons 
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(e) (f) 

FIG. 5: A simulation of the ID variable distributions (minus their requirement value) for photons in a GMSB model with 
m~o = 110 GeV/c 2 . The solid line is for prompt photons, simulated as decay photons from x?'s with a lifetime of ns and 
the dashed line is for photons from long-lived x?'s with a lifetime of 10 ns. Entries to the left of the dashed vertical line pass 
the corresponding requirement. The bin at —2.8 in (d) collects the photons that have no track within the isolation cone. In 
(f) the bin at —6 shows the photons that have no 2nd CES cluster nearby. The distributions for all ID variables do not change 
significantly between the prompt and the long-lived case except for slight deviations in the energy isolation in (b) as discussed 
in the text. 
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FIG. 6: The efficiencies for photons and electrons to pass the 
ID requirements in Table [I] vs. incident angles a and /3. The 
solid squares represent MC photons from Xi ~ * lG decays 
(m~o =110 GeV/c 2 , r~o = 10 ns) while the empty circles 
represent electrons from a W — ► ev data sample that pass the 
requirements in Table[Tl] The efficiency falls by ~15% from 0° 
to 60° in p. This effect is mostly due to the energy isolation 
requirement, as discussed in the text. 



and spikes to real electrons from W — » ev events. The 
photon candidates pass all but the Ap identification re- 
quirements shown in Table [l]in events with ft? > 30 GeV, 
while the electrons selected pass the requirements in Ta- 



ble IV As shown in the figure, a requirement of Ap < 0.6 
rejects ~100% of all spikes with a minimal loss in effe- 
ciency for real photons. Thus, this source will be ne- 
glected in the background estimate. 
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FIG. 7: A comparison of the PMT asymmetry, Ap, for a 
photon+_^T sample that contains both PMT spikes and real 
photons, and a sample of electrons from W — > ev events. 
PMT spikes can be effectively removed by requiring the asym- 
metry to be less than 0.6. 



A. Measurement of the Collision Time and 
Position 

The corrected photon time is a combination of the mea- 
surements of the photon arrival time and position using 
the EMTiming system and the primary interaction posi- 
tion and time using the COT. We begin with a descrip- 
tion of a new vertexing algorithm that provides this time 
and continue with the EMTiming measurement and the 
final t corr calculation. 

The standard vertexing algorithms [22] reconstruct the 
vertex position (xi) from high quality COT and SVX 
tracks. However, it is important also to measure to and 
to separate tracks from the vertex that produced the pho- 
ton from any other vertex that lies close in space but oc- 
curs at a different time. This is particularly true at high 
instantaneous luminosities where two or more collisions 
can occur in one event and can lie close to each other in 
z. Misassigned vertex events are a dominant contribution 
to the background estimate. 

To solve this problem, we have developed a new vertex 
reconstruction algorithm based on track clustering. The 
procedure [23| uses tracks with a well measured t and z 



that pass the requirements in Table III and groups those 
that are close to each other in both space and time. The 
algorithm can be separated into three phases: (1) the 
initial assignment of tracks that are nearby in to and zq 
into clusters, (2) the determination of the to, 2rj> an d £pr 
of the vertex, and (3) the adjustment of the number of 
clusters by merging clusters that are close to each other. 

A simple algorithm is used to make a preliminary as- 
signment of all tracks into clusters. It is designed to over- 
estimate, initially, the number of vertices in the event to 
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p T > 0.3 GeV/c 

p T > 1.4 GeV/c or passes the slow proton rejection cuts 
if charge > 

M < 1-6 
\z \ < 70 cm 
Err(zo) < 1 cm 
|to | < 40 ns 

0.05 < Err(fo) < 0.8 ns 
Traverses >3 stereo and >3 axial COT 
superlayers with 5 hits each 



TABLE III: The set of requirements for tracks to be included 
in the vertex reconstruction. These are the standard track- 
ing requirements [8] , but with additional quality requirements 
on the to measurement and a slow proton rejection require- 
ment [21] to remove tracks that likely have a mis-measured 
track t . These variables are described in the Appendix. 



Er > 20 GeV and |r?| < 1.0 

Fiducial: not near the boundary, in <j> or z, of a 

calorimeter tower 
£ H ad/£Em < 0.055 + 0.00045 ■ E 

xLip < io 

Lshr < 0.2 

p T > 10 GeV/c 

E lBO < 0.1 -Er 

—3 < A.t • q < 1.5 and |Az| < 3 cm 
|zo| < 60 cm 

p T > 50 GeV/c or 0.5 < E/p < 2.0 
Track traverses >3 stereo and >3 axial COT 
superlayers with 5 hits each 



TABLE IV: The identification requirements for use in select- 
ing electrons with high purity to study the vertexing perfor- 
mance. Note that "cp is the charge of the electron. The 
identification requirements are summarized in the Appendix 
and described in more detail in [8]. 



obviate the need for dividing a single cluster into two 
separate clusters, called splitting. The highest-p T track 
is designated as the "seed" of the first cluster, and any 
lower-p T tracks that lie within three times the typical 
cluster RMS (0.6 ns and 1.0 cm for to and zo, respec- 
tively) are also assigned to it. The highest-p T track from 
the remaining set of tracks is then picked as a second seed 
and tracks are assigned to it, and so forth until no tracks 
are left. The mean position and time for each cluster, 
-^vertex and * V ertcx respectively, is then calculated. 

The vertexing algorithm is essentially a likelihood fit 
using an iterative procedure to get a best estimate of 
the true number of vertices and their parameters |24j . 
We allow the cluster parameters to float in the fit and 



maximize the probability that each track is a member of 
a vertex with a track density that is Gaussian in both 
space and time. All clusters are fitted simultaneously. If 
during the procedure the means of two clusters are within 
both 3 cm in zq and 1.8 ns in to or if two clusters share 
the same set of tracks, then the clusters are merged. No 
splitting is done because the initial seeding is designed to 
overestimate the number of clusters. Splitting a cluster 
with a too-large RMS can result in two clusters that both 
do not pass the final requirements and would reduce the 
clustering efficiency. Having two clusters merged that are 
close in both space and time does not substantially affect 
the * corr measurement. We choose the primary vertex for 
an event to be the highest Ep-r cluster that has at least 
4 tracks. 



1. Vertexing Resolution and Efficiency 

The cluster resolution, the reconstruction efficiency, 
and beam properties are measured using a high purity 
W — > ev data sample, selected using the cuts in Table |TV| 
To measure the performance for possible events with pho- 
tons, the electron track is removed from the vertexing 
and is used to measure the vertexing performance as it 
identifies the correct event vertex. Figure [8] shows the z$ 
and *o distributions as well as their correlation for the 
vertices in this sample. Both are roughly Gaussian and 
centered at zero with an RMS of 25 cm and 1.28 ns, re- 
spectively, reflecting the accelerator performance. There 
is a non-Gaussian excess around zero in the z distribu- 
tion that comes from events that contain more than one 
vertex. In this case the clustering has merged two vertices 
that are close to each other, which most likely happens 
at z = cm. The correlation between the collision po- 
sition and time distributions is caused by the differences 
in the proton and anti-proton bunch structure within the 
accelerator (a p ~ 50 cm and ap f=a 70 cm [7]). 

The vertexing resolution is estimated using a subsam- 
ple of the events with only one reconstructed vertex. For 
each event the tracks in the vertex are randomly divided 
into two groups that are then separately put through 
the vertexing algorithm. Figure [9] shows the distance 
between the two clusters, divided by v2 to take into ac- 
count the two measurements, giving a resolution mea- 
surement of at = 0.22 ns and a z — 0.24 cm. The sec- 
ondary Gaussian in Fig. (9jb) indicates cases where two 
different vertices have been combined into one cluster. 
Figure [10] shows the difference in time and position be- 
tween the reconstructed cluster and the electron track 
(not included in the vertexing) for the full sample. The 
distributions are well described by two Gaussians that are 
both symmetric and centered at zero, indicating no mea- 
surement bias. The primary Gaussian distribution con- 
tains events where the reconstructed cluster is the vertex 
that produced the electron. Its RMS is dominated by the 
resolution of the electron track position and time. The 
secondary Gaussian distribution contains events where 
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FIG. 8: Plots (a), (b), and (c) show the to, Zo, and their cor- 
relation respectively for the reconstructed highest Epr vertex 
in W — > ev events. The fits in (a) and (b) are both a single 
Gaussian. The falloff in the (b) at \z\ ~ 60 cm is due to the 
requirement that all tracks have \z\ < 70 cm. In the search 
the vertex is required to have \z\ < 60 cm. 
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FIG. 9: The difference in t and z between two arbitrarily 
selected sets of tracks from the same reconstructed vertex in 
a W — > ev dataset with the electron track removed from the 
vertexing. This is a measure of the vertex resolution, (a) is 
fit with one Gaussian while (b) is fit with two. Note that the 
factor of \/2 is already taken out. 



the electron does not originate from the highest Spy ver- 
tex in the event. 

The efficiency of the vertex reconstruction algorithm is 
investigated using two separate methods. The efficiency 
as a function of the number of tracks is determined by 
selecting events that contain a cluster with a high track 
multiplicity. Next, various random subsets of the tracks 
are taken that belong to this cluster to see if they alone 
could produce a cluster. Figure [TT] shows the ratio of 
subset samples in which a cluster is reconstructed to all 
cases tried for a given set of tracks as a function of the 
number of tracks in the various subsets. The algorithm is 
over 90% efficient if 4 tracks are present, where the ineffi- 
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FIG. 10: The difference in t (a) and in z (b) between the elec- 
tron track and the highest Epy reconstructed vertex (without 
the electron track participating in the vertexing) mW^ev 
events. The distributions are centered at zero and fit with 
double Gaussians, indicating that there is no bias in the clus- 
tering procedure. The secondary Gaussian contains events 
where the electron does not originate from the highest Spy 
vertex in the event. 



ciency is usually caused by the algorithm reconstructing 
two separate clusters each with <4 tracks, and 100% ef- 
ficient with 6 tracks (the final analysis requires at least 
4 tracks). A second method that also allows for a mea- 
surement of the efficiency as a function of the Sp^ is to 
consider tracks in a 2 cm x 2 ns window around the elec- 
tron track (~5cr in each direction) and search for clusters. 
Only events with at most one reconstructed vertex are 
considered. While this result is not biased by selecting 
cases with a known vertex, the disadvantage is that for 
resolution reasons not all tracks are in the window, re- 
sulting in a small under-counting of the number of tracks. 
Figure [TT] shows that the efficiency as a function of the 
number of tracks in the vertex yields a similar result for 



the two very different methods. This gives confidence in 
the results as a function of Sj»t- The efficiency plateaus 
at T,pT = 7 GeV/c, as higher p T tracks have a better to 
resolution measurement. It is important to note that the 
efficiency measurements as a function of the number of 
tracks are sample-dependent. For instance, if a sample is 
chosen that is biased towards a higher average track p T 
then the efficiency might be higher for a smaller number 
of tracks, or if a sample contains many high-p T tracks, the 
efficiency as function of Spy might plateau earlier. Since 
the search requires Spy > 15 GeV/c, as stated later, we 
take the efficiency for the vertex selection requirements 
to be 100%. 



B. The Corrected Photon Time 

With the vertex time and position in hand, we move 
to a full measurement of i cor r by incorporating the EM- 
Timing information. The time of arrival recorded by the 
EMTiming system TDCs is corrected using calibrations 
that take into account channel to channel variations and 
an energy-dependent ( "slewing" ) effect due to the fixed- 
threshold discriminators. A full description of the hard- 
ware as well as the correction and calibration procedure 
is described in Ref. [16J. The t COIT resolution for electrons 
from W — > ev events is 0.64 ns (0.63 ns) for collision data 
(MC), dominated by the intrinsic resolution (0.5 ns), the 
precision of the TDC output (0.29 ns) and the vertex to 
resolution (0.22 ns). A comparison of detector simulation 
to collision data for W — > ev events is shown in Fig. 12 
There are no non-gaussian tails out to ~ 5c 



III. TRIGGERS, DATASETS AND EVENT 
PRESELECTION 

The event selection is a three stage process. The stages 
are (1) an online sample is selected (during data taking), 
(2) a 7+jet + ^r "preselection sample" is selected offline, 
and (3) the event selection uses optimized final event se- 
lection requirements. The full set of requirements that 
determine the preselection sample for the search are sum- 
marized in Table [V] The optimization and final event 
requirements are described in Section |VI| 

The analysis begins by selecting events online using a 
single set of 3-level trigger requirements that require a 
photon candidate and Tfrp. The Level 1 trigger requires 
a single tower in the calorimeter with \q\ < 1.1, E-p > 
8 GeV, E Had /E EM < 0.125 GeV, and > 15 GeV. 
For a description of the ID variables, see the Appendix. 
The Level 2 trigger requires the event to have an EM 
cluster with Ep > 20 GeV and fa > 15 GeV. At Level 
3 the requirements are tightened with Ep > 25 GeV, 
En a d/E EM < 0.125, and $p >25 GeV. The data con- 
sist of events from the data-taking period from Decem- 
ber 2004, when the EMTiming system became fully func- 
tional, until November 2005. The data correspond to an 
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FIG. 11: The clustering efficiency as a function of the number 
of tracks using (a) the subset method and (b) the window 
method, and (c) as a function of the EpT of the tracks using 
the window method. Note that a cluster is required to have at 
least 4 tracks and the efficiency is 100% for EpT > 15 GeV/c 
in this search. 
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FIG. 12: A comparison between MC (solid) and collision data 
(points) for t corr for electrons from a W — > ev sample. The 
distributions are well centered around and the resolutions 
of collision data and MC fit well with a fully corrected RMS 
of 0.64 ns. 



integrated luminosity of (570±34) pb 1 . 

The sample of 7 + tfsr candidate events that pass the 
trigger requirements is processed offline where the event 
characteristics are refined to increase the signal purity 
and further reduce the backgrounds. The offline pre- 
selection requirements include photon ID and pr re- 
quirements as well as jet, vertex, and cosmic ray re- 
jection requirements. To ensure that all signal events 
would have passed the trigger with 100% efficiency each 
event is required to have fir > 30 GeV, a photon with 
Ej > 30 GeV, and pass the identification criteria shown 
in Table |U 

We require the presence of at least one jet and a high 
Ep-r vertex in each event for the preselection sample. 
This preserves the acceptance of X2X1 an d X1X1 P ro ~ 
duction while maintaining a search strategy that is as 
model-independent as possible. While the term "jet" 
typically refers to the hadronization of a high energy 
quark or gluon that is produced in the collision, at CDF 
jets are identified as clusters of energy in the calorime- 
ter. Hence, the hadronic decays of taus and/or the 
energy deposits from electrons or photons are also effi- 
ciently reconstructed as jets. Requiring at least a single 
jet with Er > 30 GeV and \r]\ < 2.1 [H] retains high 
efficiency and significantly reduces non-collision back- 
grounds which typically only produce a single photon 
candidate. As previously mentioned, each event must 
also have a good space-time vertex with at least 4 good 
tracks and a Upr of at least 15 GeV /c. This allows for 
a good t corr measurement and further helps reduce the 
non-collision backgrounds. We also require \z\ < 60 cm 
and I to I < 5 ns for tracks to be included in the vertexing 
so that both the COT tracking and the calorimeter are 
able to produce high quality measurements. 

A cosmic ray that traverses the detector can create hits 
in the muon system that are not associated with tracks 
in the COT and deposit a photon candidate nearby in 
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Photon 



Er > 30 GeV and \r/\ < 1.0 

Fiducial: not near the boundary, in <f> or z, of a 

calorimeter tower 
ShW^em < 0.125 

Energy in a A_R = 0.4 cone around the photon 

excluding the photon energy: 

E Iso < 2.0 GeV + 0.02 • (Er - 20 GeV) 
No tracks pointing at the cluster or 

one track with p T < 1.0 GeV/c + 0.005 ■ Er 
TiPt of tracks in the 

AR = 0.4 cone < 2.0 GeV/c + 0.005 • Er 
No second cluster in the shower maximum detector 



A r 



p2 cluster 
-EPMT1+-EPMT2 



< 0.6 



Jet 



Ef > 30 GeV 



jet I 



< 2.0 



Highest Spy Space-Time Vertex with: 

A trks > 4 

Ep T >15 GeV/c 
|*| < 60 cm 
|£o| < 5 ns 

Global Event Cuts 

pr > 30 GeV 

Passes cosmic ray rejection requirements 



TABLE V: The requirements used to obtain the preselection 
sample of 7 + jet + fir events. The cosmic ray rejection cut 
is described in more detail in [25]. The number of events in 
the data that pass each cut are shown in Table |VI| For more 
detail on the ID variables, see the Appendix. 



the calorimeter. An event is rejected from the preselec- 
tion sample if there are potential cosmic-ray hits in the 
muon chamber within 30 degrees in (j) of the photon that 



are not matched to any track. Table VI lists the cumu- 
lative number of events that pass each of the successive 
requirements to create our preselection sample. 



IV. BACKGROUNDS 

Backgrounds to the 7 + jet + fr signature can be 
categorized into two different classes: collision and non- 
collision events. The rate that each type of background 
contributes to the final signal time-window is estimated 
solely from collision data using control samples of events 
that pass all of the final requirements excluding timing. 
We define the "kinematic sample" as the events that pass 



Selection 


No. of Observed 




Events 


Er > 30 GeV, pr > 30 GeV, 




photon ID and fiducial requirements 


119944 


Vertex with Ep T > 15 GeV/c, >4 tracks 


19574 


>1 jet with Er > 30 GeV and \r)\ < 2.0 


13097 


Cosmics rejection 


12855 



TABLE VI: Event reduction for the preselection 7 + jet + . 
sample. For the individual requirements see Table |V| 



tributions outside the timing signal region are used to 
normalize each background, which is then extrapolated 
into the signal time region. In this section each of the 
backgrounds is described, and the signal estimation tech- 
niques are outlined. 



A. Standard Model Backgrounds: Prompt Photons 

Prompt collision events dominate the sample and pop- 
ulate the region around t COT1 - = ns. As shown later, 
it is not important for this search to distinguish further 
between the various prompt photon sources. Most events 
are from jj and jj events with one jet reconstructed as 
a photon and with If)? from the mis-measurement of the 
photon and/or jet in the calorimeter. A smaller source 
is from SM W — > ev + jets events where the electron is 
misidentified as a photon and the v leaves undetected to 
cause the fsp. In both cases these events can fall into 
the large i corr signal time window due to either Gaus- 
sian fluctuations of the timing measurement or a wrong 
collision vertex selection. The latter case dominates the 
SM background estimate and is more likely at high in- 
stantaneous luminosity when there are multiple collision 
vertices reconstructed. 

To study the t cort distribution for promptly produced 
photons, a sample of W — > ev events is selected using 
the requirements described in Table |TTJ This sample is 
used for the reasons described in Section |TTJ and has 
the additional advantage that the electron track in the 
COT allows for a determination of the correct vertex. 
To mimic closely the vertexing for events with photons, 
the electron track is dropped from the vertex clustering. 
The highest-SpT vertex is chosen as the most likely to 
have produced the EM cluster (the "photon"). Figure 13 
shows the resulting i corr distribution and has a double- 
Gaussian shape. One Gaussian comes from events where 
the vertex choice is correct and the other Gaussian comes 
from events where the vertex choice is incorrect [2S1 ■ 



the final event requirements (summarized in Section VI 
Table IX I except the timing requirement. The i corr dis- 



Figure 13 also shows these events separated into right 
and wrong- vertex subsamples, selected by requiring a 
tight match (|z track -2i vortcx | < 2 cm and |t t rack- Vertex I < 
2 ns) and anti-match between the electron track and the 
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vertex, respectively. In both cases the distributions are 
Gaussian and centered at zero. The right vertex selection 
has an RMS of 0.64 ns, reflecting the system resolution, 
and the wrong-vertex selection has an RMS of ^2.0 ns. 
The wrong- vertex time distribution can be understood by 
combining the RMS of the time distribution without the 
vertex t and z corrections (RMS= 1.6 ns) with the RMS 
of the collision to distribution (RMS= 1.28 ns as shown 
in Fig. g: RMS wrong vcrtcx = Vl.6 2 + 1.28 2 = 2.05 ns. 
The number of events in the i corr signal region (t corl > 0) 
for prompt, SM sources can thus be estimated by simple 
extrapolation from a measurement of the timing distri- 
bution for f corr < 1.2 ns. 

The systematic uncertainty on the number of prompt 
events in the signal region is dominated by the observed 
variation in the mean and RMS of the t corr distribution 
as a function of the tfsr, jet Err, and photon Err require- 
ments. To estimate the variation, we study the i corr dis- 
tribution for samples of electrons in W + jets events for 
various electron Err, jet Err and frr event requirements 
(20 < E^ lc < 40 GeV/c 2 , 25 < Ej, ct < 40 GeV/c 2 , and 
30 < pr < 50 GeV). The results are shown in Fig. [14 
The variation in the mean is up to 0.1 ns and is conser- 
vatively rounded up to 0.2 ns. Similarly, the systematic 
uncertainty on the RMS of t carr is conservatively overes- 
timated from a fit to Fig. [l4]to be 0.02 ns and is only a 
small addition. 

For wrong- vertex assignments, there is an additional 
variation in the i corr distribution as a function of photon ?/ 
due to the incorrect time-of-flight calculation. Figure 15 



shows the mean and the RMS of the < corr distribution for 
electrons from W — > ev events where the wrong-vertex is 
selected for the timing correction, as a function of tower- 
T). We take a systematic uncertainty on the mean and 
the RMS of the wrong-vertex contribution to the icon- 
distribution to be equal to the full variation. We assign 
values of 0.33 ns and 0.28 ns, respectively, to these sys- 
tematic uncertainties, the latter arising from the largest 
variations in Fig. |15| 



B. Non-Collision Backgrounds 

The fraction of non-collision backgrounds in the kine- 
matic sample that fall in the timing signal window is 
significant. To study these backgrounds, we divide 
them into two separate sources, cosmic ray muons and 
beam related backgrounds. Cosmic ray events (cos- 
mics) come from cosmic ray muons that emit photons 
via bremsstrahlung as they traverse the detector or pro- 
duce significant ionization in a large q 2 interaction with 
the EM calorimeter. Beam halo events (beam halo) are 
caused by beam particles (mostly from the more intense 
proton beam) that hit the beam pipe upstream of the 
detector and produce muons. These muons travel almost 
parallel to the proton beam direction and shower into the 
EM calorimeter to create a photon candidate; see Fig. [T6] 
In both cases the event has significant fJr that is highly 
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FIG. 13: The t COII distribution for electrons in a sample of 
W — > ev events. In plot (a) the two Gaussians correspond 
to the cases when the highest-Epr vertex is associated to 
the electron track and when it is not. These cases can be 
separated by requiring a match (b) and an anti-match (c), 
respectively, between the vertex and the electron track in both 
space and time. 
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FIG. 15: The mean and RMS of the i CO rr distribution for 
electrons from W — > ev events, where the wrong vertex is 
picked, as a function of rj. 



FIG. 14: The mean and RMS of the t COII distribution for 
electrons from various subsamples of W — > ev + jets events 
where each entry reflects a different combination of the elec- 
tron Er, jet Er and tfsr event requirements. There are slight 
shifts as the these requirements vary. While the mean of the 
distribution is close to zero, the systematic variation on the 
mean of the primary Gaussian of the prompt time distribu- 
tion is conservatively taken to be 0.2 ns in the background 
estimates. 



correlated with the photon Er and is uncorrelated with 
any collision that might occur coincident ally at high lu- 
minosity. As cosmic ray muons interact with the detector 
and produce a photon randomly in time, their time dis- 
tribution is roughly constant over the entire calorimeter 
energy integration window range of 132 ns. Beam halo 
"photons" typically arrive a few ns earlier than prompt 
photons for geometric reasons as shown in Fig. [16| How- 
ever in this case, while the rate is lower, the photon can- 



didate can also have a i cor r of ~18 ns (and multiples 
later and earlier) if the muon was created in one of the 
bunched beam interactions that can occur every 18 ns in 
the accelerator. 

The rate at which both non-collision backgrounds pop- 
ulate the signal region is estimated from collision data us- 
ing events with no identified collision. The non-collision 
sample consists of events with a photon that passes the 
photon ID criteria listed in Table |TJ frp > 30 GeV, and 
no reconstructed vertex. This sample is used to make 
timing distribution templates from pure samples of each 
type of non-collision background. Beam halo events are 
identified by the energy deposition of the muon as it 
passes through the high rj towers of the plug hadronic 
calorimeter (|?y| > 1.1) and the central EM calorimeter 
(\v\ ^ 1-1) towers at the same </> as the photon candi- 
date; see Fig. [16] The muon deposits a small amount of 
energy in most towers along its path. Hence, we count 
the number of towers in the hadronic calorimeter with 
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FIG. 16: Illustrations of a beam halo event interacting with 
the detector. In both figures the muon path is indicated 
with an arrow, (a) A comparison of the time distributions of 
prompt collision events with beam halo photon candidates for 
three example towers in the calorimeter shows that the mean 
time changes as a function of tower rj and is always less than 
zero. The y-axes are in arbitrary units, (b) An illustration of 
how the beam halo interacts with the calorimeter. The muon 
travels through multiple towers in the hadronic calorimeter 
at high rj before hitting the electromagnetic calorimeter. 



FIG. 17: The variables used to separate cosmic and beam 
halo backgrounds in the 7 + fr sample without a vertex. 
Beam halo muons deposit energy in many HAD towers as they 
interact with the detector at high rj and many EM towers as 
they traverse the central portion of the calorimeter along the 
beam halo direction. 



tion for use in extrapolating the number of non-collision 
events in the signal time window from the control re- 
gions, we convolute the distributions in Fig. [18] with the 
RMS of the interaction time of 1.3 ns as the collision time 
is uncorrelated. As will be seen, the uncertainty on the 
rate of the number of events in the signal time region is 
dominated by the statistical uncertainty on the number 
of non-collision events in the control regions. We note 
that because of the accelerator geometry there are ~40 
times more beam halo events that occur around the re- 
gion <fi ~ 0° as can be seen in Fig. [19] This explains 
the — 15° < 4> < 15° separation requirement and will be 
further used in the final background estimate procedure. 



> 0.1 GeV and \rj\ > 1.1 (nHADTowers) and the num- 
ber of towers in the EM calorimeter with > 0.1 GeV and 
\v\ < 1-1 (nEMTowers). The results are shown in Fig. 17 
for the full non-collision sample. Cosmic ray candidates 
are easily separated from beam halo candidates. This is 
because cosmics do not deposit energy in the hadronic 
calorimeter with \r)\ > 1.1 and typically only deposit 
significant energy in a single EM tower. An event is 
identified as a cosmic if it has nHADTowers = and 
nEMTowers < 5. (Note that we also ignore all photon 
candidates with —15° < 4> < 15° as beam halo domi- 
nates there.) Conversely, beam halo events are identified 
if they have no muon stubs and have both nHADTow- 
ers > 1 and nEMTowers > 4. The t COIX distribution for 
each is shown in Fig. [18] for the entire calorimeter en- 
ergy integration window and indicates that the real col- 
lision contamination is negligible. As these events lack 
a vertex, the photon arrival time is corrected assuming 
zq = and to — in Eq. [I] To create the t cor r distribu- 



C. Background Estimation Methods 

The number of background events in the signal region 
is estimated from collision data by fitting a set of control 
regions with background timing shapes and extrapolat- 
ing into the signal time window. The i corr distribution 
shape "templates" for each background source are given 
in Figs. 13 and 18 Since a sample is defined by kine- 



matic cuts alone we can estimate the number of back- 
ground events in any potential signal time window using 
sensibly chosen control regions. Thus, we can predict 
the background rate for a large variety of final kinematic 
and timing cuts and use these estimates as part of our 
optimization procedure. 

The background prediction for the signal timing region 
for each subsamplc of 7 + jet + frp events after the kine- 
matic sample requirements is done as a two-step process 
with multiple control regions. There are a number of rea- 
sons for this. Multiple control regions arc used to get a 
robust estimate of each of the background event contri- 
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FIG. 19: The number of beam halo photon candidates as 
function of detector (j). Most photons arrive at <f> £s 0. 
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FIG. 18: The t CO rr distributions for the cosmic ray (a) and 
beam halo (b) backgrounds in the 7 + tfrr sample without a 
collision. 



butions that are hard to separate; for example the tim- 
ing region {—15,0} ns is populated by both the wrong- 
vertex backgrounds and beam-halo backgrounds. Sec- 
ondly, for many of the potential kinematics-only sam- 
ples, low statistics can bias the fit results. We define 
a set of control regions chosen such that each is largely 
dominated by a single background source, and use an it- 
erative fitting procedure to ensure that each background 
is well estimated for each kinematic requirement choice 
during optimization. 

The control regions are designed to allow for a good 
estimation of each background separately. Since the cos- 
mics rate is essentially constant in time, the time control 



region is defined to be {25,90} ns and is chosen such 
that (a) it is well above the beam halo secondary peak 
at ~18 ns and (b) it does not include the region close 
to the end of the calorimeter energy integration window 
where the event rate falls sharply. The beam halo con- 
trol region is defined to be {—20, —6} ns and is chosen 
such that (a) it contains most of the beam halo events 
but (b) stays well away from the region dominated by the 
prompt photon production. The standard model control 
region is defined to be {—10,1.2} ns. An additional re- 
quirement on this region is that the photon must have 
I <f>\ > 15°. This allows us (a) to include as much of the 
collision data as possible to get good precision on the 
ratio of right to wrong- vertex events, (b) allows for a po- 
tential signal region above 1.2 ns, and (c) removes most 
of the beam halo contamination. We note that while the 
(j) restriction is useful for estimating backgrounds, it is 
not an effective tool in improving the sensitivity. While 
the upper time limit of the signal region at 10 ns is not 
quantitatively motivated, it contains most of a long-lived 
signal on the order of nanosecond lifetimes as the time 
distribution falls exponentially (Fig. [2|. 

The background prediction for the signal timing re- 
gion for each subsample of 7 + jet + fr events after the 
kinematic sample requirements is done as a two-step pro- 
cess. In step 1, the wrong-vertex fraction and the over- 
all prompt photon rate are measured using the collision 
data control sample after a correction for non-collision 
contamination. The non-collision contamination to the 
control sample is estimated by fitting the beam halo tem- 
plate and cosmics, simultaneously in the control regions 
{—20,-6} ns (beam halo dominated) and {25,90} ns 
(cosmics dominated). The fit is extrapolated to the colli- 
sion control region {—10, 1.2} ns, where the non-collision 
contamination is subtracted off. The remaining data are 
then fit using the two single Gaussian functions as shown 
in Fig. [13] While the mean and RMS of both functions 
are fixed, the normalizations are allowed to float. After 
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fitting, the final normalization is scaled by a factor of 
12/11 to account for the removed data in the 30° slice 
around <fr = 0. The statistical error on the prediction in 
the signal region is determined by the fit. The uncer- 
tainty on the number of events in the signal time window 
is estimated by varying the collision background frac- 
tions, means, and RMS's according to their systematic 
and statistical uncertainties. We note that if we choose 
the kinematic sample to be the dataset that passes the 
preselection sample criteria (see Table Iv} , the resulting 
fraction of wrong- vertex events is (3±1)%. 

In step 2, the rate of the non-collision backgrounds 
is estimated using the entire <fi region. In the full data 
sample a simultaneous fit is performed for the normal- 
ization of the beam halo and cosmic ray backgrounds 
using the control regions, after subtracting off the ex- 
pected contamination from collision sources in both re- 
gions obtained in step 1. The uncertainties on this esti- 
mate are dominated by the statistical error on the num- 
ber of events in the control regions and the uncertainty 
on the extrapolation from the prompt background. With 
this technique, the background estimation for all sources 
is robust enough to be applied to a variety of kinematic 
requirements. This feature will be used along with the 
simulated acceptance of GMSB events for the optimiza- 
tion. 



V. ACCEPTANCES FOR GMSB EVENTS AND 
THEIR SYSTEMATIC UNCERTAINTIES 

We use MC techniques to estimate the acceptance and 
overall sensitivity to GMSB models. The sparticle prop- 
erties (mass, branching fractions etc.) are calculated with 
ISASUGRA [27] . Samples of events of GMSB processes are 
simulated according to their production cross sections us- 
ing pythia [13], a full detector simulation, as well as par- 
ton distribution functions (PDFs) [25J . All sparticle pro- 
duction mechanisms are simulated as this maximizes the 
sensitivity to the model [29 . To map out the sensitivity 
for GMSB models as a function of Xi mass and lifetime, 
MC samples are generated for 65 < m^o < 150 GeV/c 2 
and < r 7 o < 40 ns. As ^5% of the simulated events 
pass all the selection requirements, the size of the MC 
samples is chosen to be 120,000 events so that their sta- 
tistical uncertainty is ~1% and negligible compared to 
the combined systematic uncertainty. 

The total event acceptance is 

A - e= [A - e)si g nai mc x Cmc, (3) 

where the MC program is used to estimate A, the fraction 
of events that pass the kinematic sample requirements 
and to estimate e, the fraction of these events that remain 
after the i corr requirement. Cmc is a correction factor for 
efficiency loss due to the cosmic ray rejection requirement 
and is not simulated. Table IVTT1 shows the breakdown of 
the number of MC events after each of the preselection 
sample requirements in Table [V] for an example GMSB 



Requirement Events passed (A ■ e) signal mc 

(%) 

Sample events 120000 100.00 

Central photon with 



Er > 30 GeV, and 



fir > 30 GeV 


64303 


53.6 


Photon fiducial and ID cuts 


46730 


38.9 


Good vertex 


37077 


30.9 


>1 jet with Er > 30 GeV 






and \r)\ < 2.0 


28693 


23.9 


Cosmic ray rejection (xCmc) 


N/A 


23.5 



TABLE VII: Summary of the MC event reduction for a GMSB 
example point at m~o = 100 GeV/c 2 and r~o = 5 ns as a func- 
tion of the preselection sample cuts of Table [V] Note that the 
efficiency loss caused by the cosmic ray rejection requirement 
is implemented as an MC correction factor, Cmc- 



point at m^o = 100 GeV/c 2 and r^o = 5 ns, near the 
expected sensitivity limit. 

The loss of signal events due to the cosmic ray rejection 
requirement is chiefly caused by real cosmic rays overlap- 
ping the signal events and causing the requirement to fail. 
This efficiency is estimated simply to be equal to the effi- 
ciency of the requirement as measured from the preselec- 
tion sample but additionally requiring the photons to be 
within |i COIT | < 10 ns to select collision events with high 
purity. There are 12,583 events in this sample. 12,360 
events remain after the cosmic ray rejection requirement, 
giving an efficiency of Cmc = ^'583 = (^8 ± 1)%, with 
the error conservatively overestimated. 

The systematic uncertainty that enters the limit calcu- 
lation (and thus a proper optimization) is dominated by 
the potential shift of the i corr measurement for the kine- 
matic sample requirements. This along with the remain- 
ing systematic effects on the acceptance, luminosity, and 
production cross section are summarized in Table |VIII| 
The uncertainty is evaluated at m^o = 95 GeV/c 2 and 
t^o = 10 ns. The effect of varying m^o and r^o is negli- 
gible when compared to the other systematic effects. We 
next describe the estimation of these important effects. 

• Time measurement: There is an uncertainty on the 
acceptance due to the systematic variations in the 
i CO rr measurement shown in Fig. |13| Three types of 
uncertainties are considered simultaneously: (1) a 
shift in the mean of i COIT measurement, (2) a change 
in the RMS variation of the i corr measurement, and 
(3) a change in the fraction of events that have 
an incorrectly chosen vertex. The variation of the 
mean of the right (wrong) vertex i corr measurement 
has been conservatively overestimated to be 0.2 ns 
(0.33 ns) and can shift events into and out of the 
signal region. The fractional variation in accep- 
tance due to this effect is estimated to be 6.7%. 
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The fractional change in acceptance due to chang- 
ing the RMS of the i corr measurement is estimated 
to be 0.03%. The variation due to fluctuations in 
the number of additional vertices, is ~1.5% [3D]. 
Taken in quadrature the total uncertainty is 6.7% 
and forms the dominant contribution to the sys- 
tematic uncertainty on the acceptance. 

• Photon ID efficiency: As described in Section [TTJ 
the systematic uncertainty on the photon ID effi- 
ciency is estimated to be 5%. 

• Jet energy: As the event selection requires a jet 
with Et > 30 GeV a systematically mismeasured 
jet can contribute to the acceptance uncertainty. 
We use the standard CDF procedure [15] of vary- 
ing the jet energy by ±lcr of the estimated energy 
systematic uncertainty and find the resulting vari- 
ation in the acceptance to be 1.0%. 

• Initial and final state radiation: The uncertainty 
in the MC simulation of ISR and FSR effects can 
cause the photon, the jet, or the fr to be systemat- 
ically more likely to pass or fail the kinematic sam- 
ple requirements and affect the acceptance. This 
is estimated using the standard CDF procedure 
of varying the ISR/FSR parameters as described 
in [25]. The systematic variation in the acceptance 
is estimated to be 2.5%. 

• Parton distribution functions (PDFs): The produc- 
tion cross section and the acceptance have uncer- 
tainties due to uncertainty in the PDFs. The uncer- 
tainty is estimated using the standard CDF proce- 
dure of varying the PDFs within the uncertainties 
provided by CTEQ-6M as described in [28]. We find 
a relative uncertainty of 0.7% on the acceptance 
and 5.9% on the cross section. 





Relative 


Factor 


Systematic 




Uncertainty (%) 


Acceptance: 




icon measurement and vertex selection 


6.7 


Photon ID efficiency 


5.0 


Jet energy scale 


1.0 


Initial and Final State Radiation 


2..') 


Parton Distribution Functions 


0.7 


Total 


8.8 


Cross section: 




Parton Distribution Functions 


5.9 


Renormalization scale 


2.4 


Total 


6.4 


Luminosity 


6.0 



TABLE VIII: Summary of the systematic uncertainties on 
the acceptance and the total production cross section. 



• Renormalization scale: There is a systematic uncer- 
tainty of the NLO production cross section which is 
estimated using the standard technique of varying 
the renormalization scale between 0.25-g 2 and 4-g 2 
using prospino2 [3_T]. The variation of the cross 
section is estimated to be 2.4%. 



VI. OPTIMIZATION AND EXPECTED 
SEARCH SENSITIVITY 

The sensitivity to sparticle production is estimated in 
the form of the expected 95% C.L. upper cross section 
limits (egg 13 ) for various points in parameter space. Be- 
fore unblinding the signal region in the data we optimize 
the search sensitivity and determine the best event selec- 
tion requirements for a prospective GMSB signal. This 
is done using the background rates and the signal ac- 
ceptances for all sparticle production, with uncertainties, 
available for different sets of selection requirements. The 
procedure is to consider the number of events "observed" 
in a pseudo-experiment, A'obs, assuming no GMSB sig- 
nal exists, and to calculate CT95(-/V bs) using a Bayesian 
method with a constant cross section prior [32 . The 
uncertainties on the signal efficiencies, backgrounds, and 
luminosity are treated as nuisance parameters with Gaus- 
sian probability distributions. We write <795(-/V bs, cuts) 
since the limit is also a function of the number of pre- 
dicted background events and A ■ e, where both factors 
depend on the set of requirements (cuts) used. 

The expected cross section limit is calculated from 
CTg5(iV bs, cuts) and takes into account the outcomes 
of the pseudo-experiments determined by their relative 
Poisson probability [33 j; V. The expected cross section 
limit and its RMS are given by: 

oo 

°95 P ( cuts ) = ^2 095 (-Nobs, CUts) 
iV obs =0 

■P(AA obs ,AA back (cuts)) (4) 

OC 

RMS 2 (cuts) = (095 (AUs, cuts) - a^ p (cuts)) 2 

JV obB =0 

• P(A/" obs) A/back (cuts)), (5) 

where Ab ac k(cuts) is the number of expected background 
for a given set of cuts and 7 5 (A b s ,A bac k(cuts)) is the 
normalized Poisson distribution of AT b s with a mean 
Ab ac k(cuts). The expected maximal sensitivity for each 
GMSB parameter choice is found when the set of require- 
ments minimizes c7g5 P (cuts). To find the minimal o~g^ p 
we simultaneously vary the photon Er, fr, and jet Err 
thresholds, A0(^r, jet), and the lower limit on t CO r r . Here 
A<fi(ftrp,jet) is the azimuthal angle between fip and the 
highest-£r jet. This angle cut helps reject events where 
the Ifsr is overestimated because of a poorly measured 
jet. The upper limit on i corr is kept constant at 10 ns. As 
an illustration of the optimization, Fig. [20] shows the ex- 
pected cross section limit for a GMSB example point [T5] 
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FIG. 20: The expected 95% C.L. cross section limit as a func- 
tion of the lower value of the t COII requirement for a GMSB 
example point with m~o = 100 GeV/c 2 and r~o = 5 ns. The 
values of the kinematic sample requirements are held at their 
optimized values. 



at m^o = 100 GeV/c 2 and r^o = 5 ns as a function of 
the lower t corr requirement. All other requirements are 
kept fixed at their optimized values. This point is close 
to the boundary of the exclusion region. 

In the region 65 < m^o < 150 GeV/c 2 , < r^o < 40 ns 
the optimal cut values have negligible variation except for 
a small variation in the optimal jet Ej< requirements and 
the lower limit on t corr . A single fixed set of final re- 
quirement values is chosen since, far from the expected 
exclusion boundaries, this results in at most a 4% loss 
of sensitivity. The final values are Photon £^>30 GeV, 
jet £r>35 GeV, A0(|^,jet)>l.O, |r>40 GeV, and 2 < 
£corr < 10 ns. For m^a = 100 GeV/c 2 and r^o = 5 ns we 
find an acceptance of (6.3±0.6)% 



IX 



gives more 



o. Table 

details on the acceptance reduction as a function of the 
requirements. Our fit to the data outside the signal re- 
gion predicts total backgrounds of 6.2±3.5 from cosmic 
rays, 6.8±4.9 from beam halo background sources, and 
the rest from the standard model with a measured wrong- 
vertex fraction of (0.5±0.2)%. Inside the signal region, 
2 < t COTT < 10 ns, we predict 1.25±0.66 events: 0.71±0.60 
from standard model, 0.46±0.26 from cosmic rays, and 
0.07±0.05 from beam halo. Table |Xl shows the various 
possible number of "observed" events and their probabil- 
ity in the no-signal hypothesis. We find for this point in 
parameter space <7g5 P = 128 fb with an RMS of 42 fb. 
The total sparticle production cross sections, <T pxo d, are 
calculated at next to leading order (NLO) by multiplying 
the LO production cross section from pythia [TB] by the 
theoretical K-factors from [33] (~1.2 for this mass range). 
A total sparticle production cross section of 162 fb is 
predicted for this point, and thus we expect to exclude 
it. A total of 5.7±0.7 signal events is expected for this 
mass/lifetime combination. 



Preselection Sample 
Requirements 

E\ > 30 GeV, f)r > 30 GeV 
Photon ID and fiducial, \rj\ < 1.0 
Good vertex, 

£t r ac ks P T > 15 GeV/c 

|r; jct | < 2.0, > 30 GeV 

Cosmic ray rejection 
Requirements after 
Optimization 

pr > 40 GeV, > 35 GeV 

A.<p(pr, jet) > 1.0 rad 

2 ns < icon < 10 ns 



Individual Cumulative 
Efficiency Efficiency 



(%) 

54 
74 

79 
77 



92 



33 



(%) 

54 
39 

31 
24 
23 



21 
18 
6 



TABLE IX: The data selection criteria and the total, cumu- 
lative event efficiency for an example GMSB model point at 
m~o = 100 GeV/c 2 and r~o = 5 ns. The listed requirement 
efficiencies are in general model-dependent. The good vertex 
requirement (95% efficient) includes the \zo\ < 60 cm cut. 
The efficiency of this cut, as well as that of the photon fidu- 
cial and cosmic ray rejection cuts, is model-independent and 
estimated from data. 



N obs a 95 {N obs ) (fb) Probability (%) 



79.9 
120 
153 
196 
239 
280 



28.7 
35.8 
22.4 
9.32 
2.91 
0.729 



TABLE X: The 95% C.L. cross section limit as a function 
of the hypothetically observed number of events, the Poisson 
probability for the number of events based on the no-signal 
hypothesis (1.3 events expected) at an example GMSB point 
of m~o =100 GeV/c 2 and r~o =5 ns, and the requirements 
listed in Table [IX] We find for this point in parameter space 
CTg5 P = 128 fb with an RMS of 42 fb. A total sparticle pro- 
duction cross section of 162 fb is predicted for this point, and 
thus on average we expect to exclude it. 



Control Region Dominant Background 


Observed Events 


— 20 < tcorr < —6 ns Beam halo 


4 


-10 < tcorr < 1.2 ns SM 


493 


25 < tcorr < 90 ns Cosmics 


4 



TABLE XI: The observed number of events in each control 
region after all the optimized kinematic sample requirements. 
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corr 



(b) 

FIG. 21: The t COIT distribution including the control and sig- 
nal regions, after all but the timing cut for all backgrounds, 
for the expected signal and the observed data. A total of 508 
events is observed in the full time window. The two observed 
events in the signal region, 2 < t COII < 10 ns, are in the first 
signal time bin (marked with a star). This is consistent with 
the background expectation of 1.3±0.7 events. 



dieted and observed events. Table IXIfl shows the results 
as each of the optimized requirements is applied sequen- 
tially along with the expectations for a GMSB example 
point. The large fractional errors on the backgrounds are 
due to the systematic uncertainty on the mean and RMS 
of the SM distributions as discussed in Section Hv] The 
large fractional errors on the beam halo and cosmic ray 
estimates are primarily due to the small number of events 
in the control regions. Neither is a problem in the final 
analysis as the absolute number of background events is 
small in the signal region. After each requirement, sparti- 
cle production would have increased the number of events 
observed in the signal region above the background lev- 
els. However, there is good agreement between the back- 
ground prediction and the number of events observed in 
all cases. The bulk of the beam halo and cosmics back- 
ground are rejected by the timing requirement. 

There are 2 events in the final signal region, 2 < 
^corr < 10 ns, consistent with the background expecta- 
tion of 1.3±0.7 events. Figure [2l|b) shows in detail the 
time window immediately around the signal region. The 
data is consistent with background expectations. The 
two events have t corr of 2.2 ns and 2.6 ns respectively. 
Figure [22] shows the distributions for the background and 
signal expectations along with the data as functions of 
the photon Er, jet Ep, tfr, and A(f>(f)r, jet) requirements. 
There is no distribution that hints at an excess. 

A model-independent exclusion limit can be assigned 
based on this non-observation. The two observed events 
and the background and its uncertainty give a 95% C.L. 
upper limit (Ngl? s ) on the number of events produced 
of iVgg s = 5.2 events. Any model of new physics that 
predicts more than this number of delayed 7 + jet + fr 
events is excluded. To make our results useful for future 
model builders to calculate cross section limits for other 
acceptance models, we calculate a correction factor, C sys , 
that takes into account the systematic uncertainties on 
the acceptance, efficiency and the luminosity, which are 
also fairly model-independent. Using the relation 



~obs 
°95 



JV 95 °sys 



C-(A- 



(6) 



and the methods to calculate <Tg}? s , we find Ng£ s ■ C sys 
5.5 events. 



VII. 



DATA, CROSS SECTION LIMITS AND 
FINAL RESULTS 



After the kinematic requirements (Table IX I 508 
events remain in the data sample. Table |XI| lists the 
number of events observed in the three control regions. 
Figure [21] shows the i corr distribution from data along 
with the signal expectations and the background shapes, 
normalized using the control regions. 

Since the number of events in the timing window 
1-2 < i C orr < 10 ns is predicted by the background es- 
timation techniques we can compare the number of pre- 



Cross Section Limits and Exclusion Regions for 
GMSB Production 



To compare our results to GMSB models we calculate 
the 95% C.L. upper limits and compare to GMSB pro- 
duction cross sections. To allow for a more detailed com- 
parison to production cross sections for any other model 
that predicts heavy, long-lived, neutral particles that pro- 
duce the 7 + jet + tfrr final state (35] we parametrize the 
acceptance using variables that are largely independent 
of the GMSB specific dynamics. 
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Requirement 




Expected Background 




Expected 


Data 




SM 


Beam Halo 


Cosmics 


Total 


Signal 




Photon, fjT, jet pre- 














selection cuts and 














1.2 < t COTI < 10 ns 


490.74±295.40 


0.27±0.12 


1.30±0.49 


492.3±295.4 


11.7±1.4 


398 


pr > 40 GeV 


162.96±76.19 


0.24±0.12 


1.17±0.46 


164.4±76.2 


10.2±1.2 


99 


Jet Er > 35 GeV 


154.52±72.96 


0.12±0.08 


0.79±0.37 


155.4±73.0 


9.4±1.1 


97 


A0($T,jet) > 1.0 


13.07±11.57 


0.10±0.07 


0.52±0.30 


13.7±11.6 


8.5±1.0 


8 


2 < fcorr < 10 nS 


0.71±0.60 


0.07±0.05 


0.46±0.26 


1.3±0.7 


5.7±0.7 


2 



TABLE XII: Summary of the expected and observed number of events from the background estimate after the preselection 
sample requirements and each requirement from the optimization, separated for each background, and the expected number of 
signal events. The expected signal numbers are for a GMSB example point at m~o = 100 GeV/c 2 and t~o = 5 ns. Note that 
the additional requirement 1.2 < t COII < 10 ns is applied at the top line to allow the background estimation methods to use the 
prompt control region to make predictions at each stage. The preselection sample cuts are listed in Table |V| The background 
predictions match well with the observed number of events for each requirement indicating the background estimation methods 
are reliable. There is no evidence of new physics. 



There are several effects that cause the acceptance to 
vary as a function of both the Xi mass and lifetime. The 
dominant ones are the probability that (a) at least one 
Xi of the two decays in the detector volume to produce 
a photon that passes the kinematic sample selection cri- 
teria (P vo i) and that (b) i corr is within the signal time 
window (Pt). We find these are roughly independent of 
each other, and define A ■ e = P vo \ ■ P t ■ P CO it, where P CO rr 
is a minor correction described below. We find: 

P vol = (-0.254 + 6.85 • nr 3 TO 9 o - 1.54 • 10" 5 m| o ) 

-0. 625 + 0. 0647-m_ 

.(l_e~ ^? + °- 842 (7) 
P t = (-0.0449 + 8.69 • 10~ 3 m ? o - 3.49 • 10" 5 m| o ) 



■(l-(l-e^ +1 - 21 ) 2 ) (8) 

where each function consists of two multiplicative terms: 
a mass-dependent term that determines the overall scale, 
and a lifetime dependent term. Here m^o is in GeV/c 2 
and t^o is in ns. The small mass dependency of the over- 
all scale and of the exponential term in P vo \ both come 
from variations in the Xi boost with its mass in produc- 
tion [2] . A higher xl boost can cause the xl to leave the 
detector with a higher probability given its lifetime and 
cause the photon to be emitted at smaller angles relative 
to the Xi direction such that its arrival time becomes sim- 
ilar to a promptly produced photon. A variation in the 
boost is caused by a change in the shape of the px distri- 
bution as a function of the Xi mass. Another important, 
but non-dominant, factor is the lifetime term in the de- 
nominator of both exponentials. This takes into account 
the effect that both the acceptance and efficiency are not 
zero at low xl lifetimes but have a finite contribution due 
to the resolution of the i corr measurement. This causes 
prompt photons to fluctuate into the signal time win- 
dow. An additional lifetime dependent correction term, 



Pcom is introduced to compensate for remaining small 
deviations in A ■ e: 



Pcorr = 1.04 - T 9 ^ — — , (9) 

55.0 Xl 0.06 + (l-r ? o) 2 ' 1 ' 



where r^o is in ns. This simple parameterization well 
characterizes the acceptance for any GMSB model to bet- 
ter than 4% and gives us confidence that it can be of use 
to future model builders. 



Figure 23 shows the expected and observed cross sec- 
tion limits along with the NLO production cross section 
as a function of Xi lifetime at a mass of 100 GeV/c 2 and 
as a function of Xi mass at a lifetime of 5 ns, close to the 
limit of the expected sensitivity. Indicated is the 6.5% 
uncertainty-band on the production cross section. The 
band also shows the ±1<t statistical variations of the ex- 



pected cross section limit. Figure 24 shows the contours 
of constant 95% C.L. cross section upper limit based on 
the two observed data events and has its best sensitivity 
for lifetimes of ~5 ns. Figure 25 shows the 95% C.L. ex- 



clusion region for er pro d > 095 and <7 P rod > Cgs" 3 - Since 
the number of observed events is above expectations, <jg§ s 
is slightly larger than Ug5 P . The Xi mass reach, based on 
the expected (observed) number of events, is 108 GeV /c 2 
(101 GeV/c 2 ) at a lifetime of 5 ns. There is no exclusion 
of GMSB models with x? lifetimes less than ~1 ns as only 
few of the Xi have a long enough lifetime to produce de- 
layed photons. However, most of the parameter space 
there is already excluded by searches in 77 + fr ITU] . 
The large mass limits extend beyond those of the LEP 
searches [TO] (using photon "pointing" methods) and are 
currently the world's best. 
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P T (GeV) A<|)(g T , Jet) 



(c) (d) 

FIG. 22: The predicted and observed photon Er, jet Er, pr, and A^($r,jet) distributions for the signal region after the final 
event selection requirements. The GMSB distributions are for m~o = 100 GeV/c 2 and r~a = 5 ns. There is no evidence for 
new physics. 



B. Future Prospects 

This search extends the exclusion region close to the 
most important region of GMSB parameter space where 
the G is predicted to be thermally produced in the early 
universe with a mass of 1-1.5 keV/c 2 [13.. With a higher 
luminosity this search technique will be sensitive to this 
mass range. To investigate the prospects of such a search 
we calculate the expected cross section limit assuming, 
for simplicity, that all backgrounds scale linearly with 
luminosity (the uncertainties remain a constant fraction 
of the background). While this assumption allows for 
a quick estimate, it does not reflect the probable im- 
provements in the background rejection methods or the 
worsening effects due to the higher instantaneous lumi- 
nosity that could cause a higher fraction of background 
events with a wrong- vertex selection. As these effects 



would tend to balance each other, it can be considered 
to provide a reasonably balanced estimate. The resulting 
cross section limit improvement, along with the expected 
95% C.L event limit, Ng^ p , are shown in Table 
our example point at m^o = 100 GeV/c 2 and t^o 



XIII 



for 



< v 1 < ■ i w w i . < hum ■, • ' — o ns . 

1 . 1 

Figure [26J shows the expected exclusion region for a lumi- 
nosity of 2 and 10 fb _1 along with the parameter space 
where 1 < < 1.5 keV/c 2 . The figure suggests that 
this search technique will be sensitive to all of this im- 
portant parameter space at 10 fb _1 luminosity for Xi 
masses of less than ~140 GeV/c 2 and lifetimes of less 
than 30 ns. 
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FIG. 23: The expected and observed cross section limits as 
a function of the lifetime at a mass of 100 GeV/c 2 (a) 
and as a function of the mass at a lifetime of 5 ns (b). 
Shaded green (darker shading) is the 6.5% uncertainty-band 
for the production cross section. The yellow shaded region 
(lighter shading) is the variation in the expected limit due to 
the statistical variation on the number of background events 
in the signal region (~30%). 



FIG. 24: The contours of constant 95% C.L. cross section 
upper limit for the observed number of events in the detector. 
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FIG. 25: The expected and observed 95% C.L. exclusion re- 
gion along with the most stringent published LEP limits from 
ALEPH [10]. The highest mass reach of 108 GeV/c 2 (ex- 
pected) and 101 GeV/c 2 (observed) is achieved at a lifetime 
of 5 ns. 



VIII. CONCLUSION 

We have presented a search for heavy, long-lived neu- 
tralinos that decay via 7G in a sample of 7 + jet + fir 
events from pp collisions at ^fs — 1.96 TeV using the 
CDF II detector. Candidate events were primarily se- 
lected based on the delayed arrival time of the photon 



at the calorimeter as measured with the newly installed 
EMTiming system. In 570 pb" 1 of data collected dur- 
ing 2004-2005 at the Fermilab Tevatron, two events were 
observed, consistent with the background estimate of 
1.3±0.7 events. As the search strategy does not rely 
on event properties specific to GMSB models, any de- 
layed 7+jet + ^r signal (that passes our kinematic sam- 
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Luminosity 
(fb- 1 ) 


Expected 
Background 


Factor of 
Improvement 

On CTexp 


/V cxp 

JV 95 


0.570 


1.3±0.7 (2) 


1 


4.6 (5.5) 


2 


4.3±2.3 


0.46 


7.4 


10 


21.9±11.6 


0.0308 


24.8 



TABLE XIII: The expected search sensitivity improvement 
for various luminosities for a GMSB example point at 



100 GeV/c 2 



and t~o 



5 ns assuming all backgrounds 



and their uncertainty fractions scale linearly with luminosity. 
The numbers in parentheses reflect the observed values in this 
search. The resulting exclusion region is shown in Fig. [26] 
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FIG. 26: The expected 95% C.L. exclusion region after a 
scaling of the background prediction and the uncertainties 
for a luminosity of 2 fb -1 and 10 fb -1 respectively. The 



shaded band shows the parameter space where 1 
1.5 keV/c 2 . 



< 



pie cuts) is excluded at 95 % C.L if it produces more 
than 5.5 events. This result allows for setting both quasi 
model-independent cross section limits and for an exclu- 
sion region of GMSB models in the Xi lifetime vs. mass 
plane, with a mass reach of 101 GeV/c 2 at r^o = 5 ns. 
These results extend the sensitivity to these models be- 
yond those from LEP II 1 0J and are the world's best at 
masses > 90 GeV/c 2 . By the end of Run II, an integrated 
luminosity on the order of 10 fb -1 might be collected, 
for which we estimate a mass reach of ~ 140 GeV/c 2 
at a lifetime of 5 ns by scaling the expected number of 
background events. 
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APPENDIX: ID VARIABLES 



In Table |XI*V| we provide a description of the identifica- 
tion variables used in this analysis for electrons, photons, 
and tracks. 
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Photon and Electron Identification Variables 



Ficucial 

-BHad/-EEM 

£; Iso (GeV) 
Ntracks 
Ep T (GeV/c) 

p2 n cluster 

A P 
E/p 

2 

XStrip 

Lshr 

Ax • q (cm) 

z (cm) 
Az (cm) 



(GeV) 



|Aces| < 21 cm and 9 < |i?cEs| < 230 cm for the calorimeter cluster centroid 
The ratio of the energy deposited in the hadronic calorimeter behind the cluster to the energy 
in the cluster as measured in the EM calorimeter 



Energy in a cone of AR = ■sj A(f> 2 + Arj 2 — 0.4 around the object, excluding the cluster energy 

Number of tracks pointing at photon cluster 

Total pt of tracks in a cone of AT? = 0.4 around the cluster 

Energy of a second EM cluster, if any, as identified in the shower-maximum detector 
|gpMTi--EpMT2l where £ an d £ PMT2 are the two PMT energies 

Ratio of the electron energy as measured in the calorimeter to the momentum as measured 
by the COT 

A x 2 comparison of the shower-maximum profile to test beam data expectations 

A comparison of the energy deposition of the electron, in adjacent towers, to expectations 

The comparison between the extrapolated track position into the shower-maximum detector 

and the measured cluster centroid position, taking into account the track charge 
The measured z position of the electron track along the beam line 
Difference between the measured electron track zo and the measured zq of the vertex 



Good Track Selection Variables 



Slow Proton rejection The dE/dx for the track as it traverses the COT is not consistent with being from a 

slow proton. We require the timing variable to be less than 20 ns. For a description 
of the relationship between the timing measurement and dE/dx see [21] . 
The measured z position of the track along the beam line 
Uncertainty on the zo measurement 
The measured time of the track origin 
Uncertainty on the to measurement 



zo (cm) 
Err(zo) (cm) 
to (ns) 
Err(to) (ns) 



TABLE XIV: A description of the identification variables used in this analysis for electrons, photons, and tracks. 
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